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Abstract 

The flow equations of the Functional Renormalization Group are applied to the 
0(A r )-symmetric scalar theory, for N = 1 and N = 4, in four Euclidean dimensions, 
d = 4, to determine the effective potential and the renormalization function 
of the field in the broken phase. In our numerical analysis, the infrared limit, 
corresponding to the vanishing of the running momentum scale in the equations, 
is approached to obtain the physical values of the parameters by extrapolation. In 
the N = 4 theory a non-perturbatively large value of the physical renormalization 
of the longitudinal component of the field is observed. The dependence of the field 
renormalization on the UV cut-off and on the bare coupling is also investigated. 

I. INTRODUCTION 

An important technique in continuum field theory, is the functional renormalization group 
which represents a powerful method to approach both perturbative and non-perturbative 
phenomena; it is based on the infinitesimal integration of momentum modes from a path in- 
tegral representation of the theory with the help of a Wilsonian momentum cutoff [3|. The 
resulting functional flow equations interpolate between the microscopic theory at short dis- 
tances and the full quantum effective theory at large distances. 

In the past years various realizations of the functional renormalization group were devel- 
oped and, among the most renowned, there are Polchinski's and Wetterich's formulations (for 
reviews see Q j^]). At the same time, the derivative expansion, a powerful approximation 
scheme that relies on a small anomalous dimension of the field, was introduced to reduce the 
full flow equation to a small set of treatable partial differential equations jsl, 0], More recently, 
a different approximation scheme was presented in js-lO], where the flow, rather than being 



projected on the semi-local derivative expansion of the action, is projected on the n-point 
1-particle irreducible (1PI) vertices, so that the momentum dependence of the correlation 
functions is properly taken into account. 

The range of application of the functional renormalization group is wide and it is especially 
suitable for the study of phase transitions, due to its flexibility even in the presence of strong 
correlations or couplings. For instance fixed point studies of Ising-like or O(N) theories, 
including the determination of universal critical indices, were repeatedly carried out, using 
different versions of the renormalization group (RG) flow equation (see for instance and 



references therein), both to analyze the properties of the theory and to test the accuracy of 



the method through a comparison with other field-theoretical techniques (see e.g. 12j-[l4|) 



Some attention has also been devoted to the study of the non-perturbative features of the 
effective action and effective potential in the presence of spontaneous symmetry breaking. 
In fact is it a well known result of quantum field theory 15Ml7| that the effective potential 
defined through a Legendre transformation must be convex, a property that, in the limit of 
infinite volume, leads to a non-analytic behavior with a typical degenerate flat region of the 
potential in correspondence of the classically forbidden region between the classical minima. 
The convexity of the potential can only be recovered by resorting to some non-perturbative 
scheme (see e.g. [181 ]) and it has already been studied by means of RG techniques 0-3. 

In this paper we analyze in detail the spontaneously broken phase of the 0(iV)-symmetric 
scalar theory, for N = 1 and N = 4, in four Euclidean dimensions, d = 4, going beyond the 
Local Potential Approximation (LPA) and including the effects of the field renormalization. 
Instead of deriving the flow equations for the potential V and the field renormalization Z 
directly from an expansion of the full flow equation for the effective action around a constant 
field configuration, we shall derive them starting from the approach developed in js iO] and 
adopting the regulator for the infrared modes introduced in 2jj-[32j, which has the specific 
property of optimizing the determination of the critical exponents in the LPA. This particular 
choice allows us to analytically perform the integration over the momentum variable, which 
simplifies the structure of coupled Partial Differential Equations (PDEs) for V and Z. The aim 
of the present analysis is to provide a better understanding of the details of this phenomenon 
which, for N = 4 and d = 4, could provide some indications on the description of the broken 
chiral symmetry phase for a two-flavor quark model or on the mass generation mechanism in 
the Standard Model, although in both cases the similarity has serious limitations. In fact, in 
the former case the 0(4) model is in the correct universality class to describe the universal 
features of the chiral phase transition, but when the system is not close to the critical line, 
which is the case considered here, then some properties of the fermionic model could be 
poorly reproduced by the scalar fields. Concerning the latter case, the coupling to gauge fields 
is turned off in the 0(N) model and therefore the massless Golstone bosons play an important 
role in the infrared sector of the ordered phase and this could lead to a different behavior with 
respect to the usual Higgs mechanism where the Goldstone are replaced by the longitudinal 
degrees of freedom of the massive gauge bosons which, due to their mass, are less relevant for 
the infrared dynamics. 

An attempt to use the two coupled flow equations for V and Z in this context was already 
done in 24J for the simple quanto-mechanical problem (d = + 1 dimensions) and in 27] for 
the N = 1 scalar theory in d = 4 dimensions. In both cases the Proper Time version of the 
RG flow 33M38 | (which can be regarded as an approximation to the class of exact background 



field flows 2^39-411) was adopted in the exponential form, which is derived by using a sharp 
cut-off on the Proper Time variable 36|, |42( . In the quanto-mechanical case the corrections 
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induced by the inclusion of Z turned out to be in very good agreement with the (exact) 
Schroedinger equation output while, for the quantum field theory, a very large Z was observed 
in the classically forbidden region. However, in the latter case, only partial conclusions could 
be drawn, because the numerical instabilities in the flow equations made rather difficult the 
approach to the infrared (IR) limit, where all quantum modes are integrated out and the full 
field renormalization is obtained. On the contrary, the approach used in this paper, based on 
a different approximation scheme, turns out to be numerically more stable, thus allowing us 
to push the RG scale low enough to get a better insight of the IR region. 

The scheme of the paper is the following. In Sect. II the main details of the RG flow 
equations and of the approximation used are recalled; the results of the numerical analysis 
of the potential and the field renormalization for the single field, N — 1, and for the N = 4 
theory in d = 4 are shown in Sect. III. The conclusions are reported in Sect. IV. 



II. RG FLOW EQUATIONS 

The starting point of our analysis is the partition function of the theory, defined as a 
functional of the source J(x), and also dependent on the momentum scale k: 

Z k [J] = J Vip e - s ~ ASk+ ^ J * . (1) 

where the usual action S is modified by the k dependent regulator, quadratic in the field ip, 
&S k [<p]: 

*S k [tp] = ijT R k (q)<p(q)<p(-q), (2) 

and 1=1 d d x . (3) 



with 

d d q 



<i 



(27r) d 

In order to obtain a physically relevant flow, R k (q) must suppress the modes with q ^ k 
while allowing to integrate those with q ^> k and therefore one can choose R k {q) of order k 2 
in the former region and R k (q) ~ in the latter. In particular at k = the regulator must 
vanish so that Z k=0 [J] coincides with the standard partition function Z [J]. In this framework 
the following flow equation is obtained \^ (d t = kd k ): 



1 



</ 



dtTM = - I d t R k (q) r?[q, -q; <P] + R k (q) (4) 



where T k [(p] is the modified effective action at scale k, defined by: 

/ J0 - AS k [0] - T k [0] = log Z k [J] (5) 

J X 
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(2) 

with 4>(x) = 8\nZ k [J\/8J(x) and and T k '[q,—q;(p\ is the Fourier transform of the second 
functional derivative of [</>]: 



Note that the term 



G k (q,<P)= r^[ q ,-q;<P]+R k (q)Y [ (7) 



appearing in Eq. (jlj) is the full propagator of the action modified by the regulator AS k . Then 
one can take the initial condition for the flow equation Eq. (j4]) at the scale k = A where 
fluctuations are frozen by AS k , so that r*. = A[0] « S[<f>] and the full effective action r[0] of the 
original theory is obtained as the solution of Eq. (jlj) when k — > 0, where Rk{q) vanishes and 
all fluctuations have been integrated out. 

The most straightforward approximation scheme to treat the RG flow is to express the 
effective action T k [(f)] through an expansion in the derivatives of the field: 

r*M = J (vm) + \z k {<j>) (d<pf + o(<9 4 )) . (8) 

and to project the full equation (jlj) onto a set of flow equations for the coefficients of the 
expansion Vk(4>), Z k (4>), ... . This projection requires an expansion of the right hand side of 
Eq. (jlj) in powers of the field derivatives or, equivalently, in powers of the momentum. 

An alternative scheme has been developed in j8-10] and, below, we shall briefly recall its 
essential features (in the rest of this Section we will follow the notation adopted in [lOj ). This 
scheme consists in an expansion in terms of the n-point functions, T k n \ i.e. the n-th functional 
derivative of r k [4>} with respect to the field <j). For instance the first two equations for n = 
and n — 2, correspond to the two following equations, respectively for the potential V k defined 
in Eq. (jHj) and equal, up to a volume factor, to T k evaluated at a constant field configuration 
<fi, and the 2-point function computed at constant (p : 

d t V k (<t>) = ~J^d t R k (q)G k (q,<t>), (9) 

d t T k 2 \p, 0) = Jd t R k {q)G k {q, 0) jrf(p, q, -p - q, <P)G k (q + p, <f>)T® (-p,p + q, -q, 0) 

-^r[ 4) (p,-p,g,-g,0)|G fc (g,0). (10) 

It is easy to check from the functional structure of Eq. (jlj) that the flow equation of 
involves the s-point functions with s < n + 2. Therefore an approximation or truncation is 
needed to obtain a closed set of PDEs. As explained in js iO] this can be realized at any fixed 



4 



n, by neglecting the internal momentum q in rj^ +1 ) and T k n+2 \ so that they can be written 
in terms of rj^ according to the following property, valid for constant field configurations : 

rl n+1) (te},o,0) = ^rf(fe},0) (ii) 

where the index i runs between 1 and n. In particular, for n = 2, F k and in Eq. ( flOl) are 
computed at q = and expressed in terms of T k , according to Eq. ( TTTj) . Thus, a closed set 
of two equations for and if is explicitly obtained in j8-10|. 



(2) 

Our aim is to extract the coupled flow equations of V*. and Zf. from those of and T, 



k 



This is easily obtained by performing a derivative expansion of the latter, to order 0{p ). To 
this purpose we restrict the form of 1^ according to the ansatze in Eq. (jHJ), i.e.: 

if V 0) = z^y + to) + o( P 4 ) (12) 

where the prime indicates the derivative with respect to the field 0. Then, the flow equation 
of Vk is given in Eq. ([9]), while the insertion of Eq. (112j) into Eq. (flOl) gives 

Z fc (0) + 0(p 4 ) = J 3 (p, 0) [Ff(0) + p 2 Z((p, 0)f - 

W) (TO)) 2 " JW) P 2 ^'(P, 0) + 0( P 4 ) , (13) 

where we have used the notation of llOl : 



J n (0) = J n (p = O,0) and J n (p,<P)= [ d t R k (q)G k (p + q,<l ) )G n k - 1 (q,<l ) ) (14) 

Jg 

with Gk given in Eq. ([7]) and 1^ in Eq. (112j) . Finally, to obtain the flow equation of Zfc 
one has to expand the integral J3Q0, 0) in powers of p 2 and systematically neglect the 0(p 4 ) 
terms in Eq. (TlBl and, as the terms proportional to p° in the right hand side vanish, one reads 
the equation for Zk(4>) from the 0(p 2 ) terms. Then, the PDEs coming from the derivative 
expansion to order 0(p 2 ) of the Tf equation derived in jsHIoj, contain the same approximation 
made on the latter, namely the neglect of the internal momentum q dependence of if, if . 

The extension to the 0(A r )-symmetric scalar theory is straightforward and it is also illus- 



trated in lOj. According to the symmetry of the theory the propagator can be written in 



terms of its longitudinal (L) and transverse (T) components with respect to the external field 



where 



G ab (p 2 A) = G T (p 2 ,p) ( 5 ab - + G L (p 2 lP )^. (15) 



P = -<f>a4>a (16) 

If the 2-point function is parametrized as 

r i?(P ; ~P, P) = T A(P, p)S a b + 0a06r B (p, p) (17) 
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and and T B are expressed in the following way (here the dots indicate derivatives with 
respect to p and the script k is omitted for simplicity) 

T A (p,p) = Z A (p)p 2 + V and T B {p, p) = Z B (p)p 2 + V , (18) 

then the following relations hold 

G T \p,p) = T A (p lP ) + R k (p) = Z T (p)p 2 + V + R k (p) (19) 

Gl\p, p) = T A (p, p) + 2pT B {p, p) + R k (p) = Z L (p)p 2 + V + 2pV + R k (p) (20) 
with Zt and Zl defined as 

Zt = Z A and Zl = Z A + 2pZ B • (21) 

The flow equations for V, Z A and Z B are 

dtV(p) = \{(N- ^(P) + Ii L (p)} (22) 

p 2 d t z A ( P ) + o{p') = 2 P {j^{p 2 z A + v) 2 + 4 l (p 2 z b + V) 2 - (l£ T + IJ L )V 2 } 

-\l% L p\Z A + 2pZ A ) - \ll T p 2 ((TV - 1)Z A + 2Z B ) + 0(p 4 ) , (23) 

p 2 «9 t Z B (p) + 0(p A ) = JP(N - l)(p 2 Z B + V) 2 - J 3 LT (p 2 Z A + V) 2 - jJ l (p 2 z b + V) 2 
+4 L {(p 2 Z A + 2p 2 Z B + 3V) 2 + 4p(p 2 Z B + d p V) (p 2 Z A + 2p 2 Z B + 3V) 

+4p 2 (p 2 Z B + d p V) 2 ] - \l? T {N - l)p 2 Z B - \l^ L p 2 (hZ B + 2pZ B ) 
- ((N - 1)/ 3 TT - J 3 LT - J 3 TL ) t> 2 - li L (3V + 2p6> p t>) 2 + p 2 Z B (/ 3 LT + J 3 TL ) + 0(p 4 ) (24) 

with the following definitions of the integrals (n > 1 and a, (3 stand either for L or T) : 

J?(p,p)= [d t R k (q)G n a - 1 (q, P )G^p + q ,p) and = (p = 0, p), (25) 

Jq 

The only undefined ingredient in the above flow equations is the explicit form of the regula- 
tor. There are various functional forms of R k (p), tested in many specific problems and whose 
properties have been studied in detail. From a practical point of view, since we are interested 
in the numerical resolution of a set of PDEs, it is preferable to take a particular regulator 
which allows us to solve analytically the momentum integrals. The cutoff 29h32| : 

Rk{q) = (k 2 - q 2 ) e(k 2 - q 2 ) (26) 
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(where indicates the Heaviside step function) makes the resolution of the integrals particu- 
larly simple and therefore we use this regulator in our analysis. 

Regarding this choice, some comments are in order. In fact the flow equations contain 
the derivative of the regulator with respect to the scale k and, as discussed above, the flow 
of Zfc is obtained after an expansion in powers of the external momentum p which, again, 
introduces derivatives of the regulator. The regulator in Eq. f[2T)j) has the form xQ(x) so that 
its derivative dtRk{q) generates two terms : the first one, (dtx) G(x), which produces relevant 
contributions and the second one x(d t Q(x)) which produces a term proportional to xS(x) 
that vanishes under integration over the momentum q, as long as no pathologies appear in the 
integrand. Therefore only the first term is to be retained in the flow equations. 

Then, when deriving the PDE for Zk from the flow of the 2-point function, an expansion in 
the momentum p of the integral Jz(jp, 4>), defined in Eq. 014p . to order 0(p 2 ) is necessary and 
this generates new terms proportional to B and 5-functions. Also, a derivative of delta-function 
of the form d Pv {y S(y)) (withy = k 2 — (p+q) 2 ) is generated, which, due to the particular form of 
the variable y, can be replaced by d qi/ (yS(y)) and the corresponding term can be calculated by 
means of an integration by parts, by recalling that the structure yS(y) gives zero contribution 
even when evaluated at the boundaries of the integral over the momentum variable q. After 
the integration by parts, by making explicit the dependence on the distributions, one finds 
(below, Ipvip) and h$l(p,q), i — 1,2,3 generically indicate regular functions appearing in the 
expansion and, again, y = k 2 — (p + q) 2 ) : 



and Eq. (1271) . evaluated at p — 0, contributes to the flow equation of Z^. The terms in 
Eq. (1271) proportional to Q(y) or Q 2 (y) give finite contributions while, by regarding the delta- 
function as the limit of a sequence of normal distributions, we discarded all terms proportional 
to yS(y) or (y5(y)) 2 or y S(y) Q(y), and calculated the last term in Eq. (l2~7j) . where the 
argument of the delta-function coincides with the upper boundary of the integral, according 



In Sect. IIHI we study the RG flow equation of a scalar theory in the broken phase in d — 4, 
concentrating on two specific cases: the N = 1 and the N = 4 theory. We use the ansatze in 
Eq. (jHJ) with the following initial condition at the ultraviolet (UV) scale k = A 



where a = 1 for N — 1, and runs from 1 to 4 for N = 4. According to Eq. 02 8p . at k = A the 
field renormalization is 1 (for the N = 4 theory Zl = Zt = 1 which, in turn means Za = 1, 



d P AJ 3 (pA) = l, u (p) + f dq 2 {h$(p,q) [Q(y)+y6(y)]} 
+ jf d( l 2 {h$(P,<l) [6(y) + y 5(y)] [Q(y) + y8(y)] + h$(p,q) 5(y)} (27) 



toJ Q Xo f(x)8(x-x ) = f(x )/2. 




(28) 
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Z% = 0) and the potential V\ only contains quadratic and quartic terms in the field, i.e. the 
only renormalizable terms in d = 4. In order to restrict ourselves to the broken phase, the 
bare mass, related to the quadratic term of the potential in Eq. ( l28|) . must be negative. In 
particular, its explicit value is fixed at —1 so that, in the following, all dimensionful variables 
are automatically expressed in units of the bare mass. It must also be recalled that the phase 
boundary is a function of the UV cutoff A and of the bare coupling A, and both parameters 
should not be increased too much to avoid a transition to the disordered phase. 



III. FIELD RENORMALIZATION IN THE BROKEN PHASE 

For practical convenience, in the following numerical analysis, instead of concentrating on 
the flow of Zk(4>) and Vk((j>), we focus on Zk{(f>) and V^'(0), i.e. on the derivative of the potential 
with respect to 0, which means that for the O(N) theory, instead of using p, we express all 
variables in terms of 



= J^J~ a = ^Tp . (29) 

and, in the following we indicate with the field at the minimum of the potential in the IR 
limit k — 0: 

= • (3°) 

The Green functions of the theory are obtained from the derivatives of the effective action 
evaluated at the minimum configuration <fi = <p and therefore <fi is crucial in the determination 
of physical observables. 

So, in N = 1 we solve the two coupled PDEs given by Eq. ( fl~3~I) and the field derivative of 
Eq. (Q, while for the N = 4 theory we solve the set of three PDEs for Z^, Zg in Eqs. fl23| [241) 
and the field derivative of Eq. ff22l) and in all cases d — 4. The numerical solution of our set 
of PDEs is obtained with the help of the NAG routine 43| that integrates a system of non- 



linear parabolic partial differential equations in the (x, t) two-dimensional plane. The spatial 
discretization is performed using a Chebyshev C° collocation method, and the method of lines 
is employed to reduce the problem to a system of ordinary differential equations. Typically 
we used Chebyshev polynomials of order 3 or 4 which already provide a stable solution. The 
routine contains two main parameters - the number of mesh points on the 'space' axis and the 
local accuracy A in the 'time' integration - that can be adjusted to control the stability of the 
solution. The first parameter is taken in such a way that the distance between two subsequent 
mesh points is Sx ~ 10 -3 , 10 -4 while A typically is taken between 10 -7 and 10 -9 . The size of 
the time step is adjusted at each integration step by the routine to keep the accuracy below 
A. The number of time steps required to converge is strongly dependent on the problem 
considered, ranging from few hundreds up to 10 7 for the hardest cases. 

The resolution of the PDEs is done with the 'space' variable constrained between < 
4> < Abound and for our purposes Abound = 4 was sufficiently large to ensure the stability 
of the solution. The required initial conditions for the PDEs at initial 'time' k = A and 



8 



0.4 0.8 1.2 1.6 

4> 



Figure 1: Vj.((f>) for N = 1 and d = 4 at various k: k = A(a), 0.7(b), 0.5(c), 0.3(d), 0.02(e), with 
A = 10 and A = 0.1. The dotted curve (f) is the one loop effective potential derivative, V/^), in 
terms of renormalized parameters. 



< cf) < (frbound are implicitly given in Eq. (128j) . The boundaries at = in the entire range 
of 'time' integration reflect the required symmetry of our model: V fc '(0) = Z' k (0) = while 
at = (fibound we keep free boundaries. Finally, the 'time' integration that starts at k = A 
should, in principle, reach k = or at least should be pushed close enough to zero that the 
solution has become invariant. But in the following analysis, due to the not smooth behavior 
of the solution, which gets worse when approaching the point k — 0, the numerical integration 
breaks down well before k = 0. Clearly this is peculiar of the broken phase because in the 
symmetric phase, where both solution for Vk and Zk are much smoother, one can easily reach 
values of the scale k where the solutions are already stable. In the broken phase the lowest 
accessible value of k depends on the bare parameters (larger values of A and A in this sense 
are preferable) and also on the accuracy A in the 'time' integration and on Sx in the 'spatial' 
mesh. However the process of improving 'space' and 'time' accuracy has the effect of increasing 
the computer time and in any case the routine seems to be unable to converge for k below 
some particular value. So the best compromise found for the minimal value of k in the present 
analysis is k = 0.02 at which the routine converges for all the examples considered. 

As a first step, we set the boundary values A = 10 and A = 0.1 in Eq. ( |28i) and analyze the 
derivative of the potential as obtained from the two coupled PDEs for N = 1 (and d = 4). The 
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results for V£(4>) are shown in Fig. [TJ Curve (a) is V^, while (b), (c), (d), (e) respectively 
correspond to V£ at k = 0.7, 0.5, 0.3, 0.02. 

Together with these RG outputs we also display in Fig. CD the derivative of the one-loop 
effective potential, (f), obtained with the same set of renormalized parameters of curve (e). 
More specifically, the one loop quantum correction to the effective potential, properly regu- 



larized by a four- momentum cut-off A, reads (see e.g. 44|) 



VM) = { a4 k+W)- Kwf k + vm) + AV " W } (31) 

where the prime, again, denotes the derivative with respect to the field and V is the bare 
potential which, in our case, is given in Eq. (128]) . In the presence of a negative square mass 
term and therefore of unstable modes, the logarithms in Eq. ( I3~i~j) develop an imaginary part 
and here we shall only consider the real part of the one loop effective potential. 

We observe that the flow equation has a different dependence on the the ultraviolet cut-off 
A with respect to Eq fl3T|) . This can be directly checked by integrating Eq. (Q in the LPA 
(i.e. fixing = 1) and also neglecting the k dependence of Vk in the right hand side of the 
equation, in order to preserve at least the leading dependence on A. In this case the flow 
equation reads: 

d k V k (<f>) = ^- 2 k2 + y„ w , (32) 

and its integration from k = A down to k = gives the following expression of the correction 
to the first derivative of the bare potential: 

V (6) V °" [A 2 2V» In (l I A2 ^ I A2V °M ) (33 ) 
to be compared with the derivative of Eq (|3ip : 

v ^=&\ 2A2 - 2V "^{ l+ vm))- (34) 

While the logarithmic term in Eqs. (|33|) and (|34|) is the same, the leading term proportional to 
A 2 has different coefficients. Eq. (13"3"|) also shows an additional subleading term. However, the 
difference between Eqs. ( 1331) and ( 1341) plays no role for the renormalized quantities because 
it can be totally compensated by a suitable choice of the counterterms so that the absence of 
divergences (i.e. the cancellation of any dependence on the cut-off A) ensures the smallness 
of the one loop corrections to the classical potential. In practice this is realized by fixing the 
proper mass and coupling counterterms in the one loop effective potential, so that its second 
and fourth derivatives computed at the minimum of the one loop potential, (pu, respectively 
coincide with the numerical values of the (right limit) of the second and fourth derivative of 
the RG generated potential at </>, as obtained from curve (e). As it is evident from Fig. [fl 
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curve (f) obtained by the above renormalization procedure, is extremely close to curve (e) 
in the region above the minimum, with the small differences mainly referable to higher loop 
corrections and to 1/A 2 terms that are neglected in our computation of the renormalized one 
loop potential. Instead, as is well known, the large difference between (e) and (f) in the 
region below the minimum is due to the failure of the perturbative expansion in recovering 
the convexity property of the effective potential which corresponds to the flat region of curve 
(e). 

Before proceeding we also quote the perturbative correction of the field renormalization as 



obtained from a one loop calculation in 16 



= TJi \ TTJTT^ t > ( 35 ) 



4tt 2 [ V o "(0) 

which has to be added to the leading term Z — 1. Eq. (I35|) . when computed at the physical 
value of the field, i.e. at the minimum of the potential 0, vanishes in the disordered symmetric 
phase (0 = 0), but it has a finite value in the broken phase, proportional to the coupling A. 
In particular, in our example with A = 0.1 , the one loop correction is Zn(<f>) = 1.69 10~ 3 . 

Let us come back to the RG generated potential. In the region of small <j> the curves b, c, d, e 



show the flattening which had been already observed many times [21|, |23M28| : the evolution 
of the potential is smooth until k reaches the infrared threshold where k 2 ~ — V^"(0 = 0) > 0, 
i.e. the region of the unstable modes which induce strong modifications to the propagator and 
therefore to the flow. When k becomes smaller than this threshold, V^(0) starts to develop a 
linear behavior in close to — (see curve b) and, when k gets smaller, the linear region 
extends to larger values of with a slope that decreases toward zero. On the other hand, 
the region with > is substantially k independent for k below threshold, with only very 
small changes when k ~ (see curves (c), (d), (e)). The region between the two regimes at 
small and large shows a sudden change of slope that becomes sharper and sharper for lower 
values of k. However when the details of this sudden change are properly enlarged, the field 
derivatives of curves (b),(c),(d) still show a continuous behavior while for the derivative 
of curve (e), obtained for k = 0.02, the change is so sharp that, according to the numerical 
precision imposed, the routine is not able to approach zero beyond k = 0.02. 

This is essentially what was expected. Indeed, it is known that at k — the potential 
must reproduce the features of the effective potential which is convex, with a flat region for 
— < < and with a non-analytic behavior, that consists in a jump of the second derivative 
of the potential at from zero to a finite value. As a consequence, the computation of the 
zero momentum Green functions from the derivatives of the effective potential is typically 
obtained by taking the limit of the derivatives from above, — > with > 0, while the 
limit from below (with < 0) is associated to tunneling processes with infinitely large time 
scales. In our numerical analysis the 'flat' part of curve (e) obtained at k = 0.02, is practically 
vanishing: it differs from zero for less than 10 -3 . Moreover, the fact that the positive branches 
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N=4 h=0A3 A=10 




0.25 0.252 0.254 

Figure 2: V£ for N = 4 and d = 4 at k = 0.02 and k = 0.01 in LPA with A = 0.13, A = 10. 



of (c), (d), (e) are practically coincident, implies that these curves provide a reasonably good 
estimate of even if the flow has not reached the point k = 0. 

We are also able to check Vk(<f)) for small and k. In fact, in this regime the relation 
V fc '(0) = —k 2 <f) has been repeatedl y d erived starting from the functional form of the flow 



equations in the LPA (see e.g. [19|, |20|, |22|, |25|, |28j). In the example displayed in Fig. [TJ the 



effect of Zk is included and we have no analytical solution for the coupled PDEs; however 
for the curves (b),(c),(d) and (e) we find : -V^<p) / (k 2 <p) = 0.955(10) at = 0.1 and 
—V k \(p)/(k 2 (j)) = 0.930(10) at = 0.5 and this gives an estimate of the size of the corrections 
to the above relation, which stay below 10% even for values of k and close to 0/2. 

In the N = 4 case the derivative of the potential shows a linear behavior in the range of 
small and k similarly to the case N = 1 in Fig. [TJ However for continuous symmetries in 
the ordered phase the presence of the Goldstone bosons has a strong impact on the infrared 
sector of the theory. In particular, according to the study of the susceptibility of the O(N)- 



symmetric theories (N > 1) [12|, |45|, |46[, in d = 4 one expects a vanishing second derivative of 
the potential at = 0. This means that the non-analyticity developed in the N = 1 potential 
at k = 0, discussed before, should disappear in N = 4. To analyze this point, we plot in Fig. 
|2] an enlarged detail of VI around 0, with N = 4, d = 4, for k = 0.02 and k = 0.01. This 
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Figure 3: Left panel: Z k for N = 1 at k = 0.5, 0.02 (dot-dashed); Z T for N = 4 at k = 0.3 
(dashed) and k = 0.02 (continuous). Right panel: Zl for N = 4 at k = 0.3 (dashed) and k = 0.02, 
(continuous), d = 4, A = 10 and A = 0.1 for all curves. 



computation is performed in the LPA at A = 0.13 and A = 10, in order to get closer to k = 
and have a more accurate check on the behavior of V£. In fact, in Fig. [2]the magnification is so 
large that the various mesh points can be identified but, both at k — 0.02 and 0.01, Vl shows 
a smooth growth without apparent jumps or rapid numerical fluctuations which were instead 
observed in the N = 1 case. According to the numerical precision considered we can state 
that the LPA solution is approaching V^ =0 (<j)) = 0. The inclusion of the field renormalization 
makes more difficult the approach to k = but there is no qualitative modification of the 
picture observed in the LPA up to k = 0.02. 

Let us now focus on the other variables, that is the field renormalization Zh and Z^, Z?- The 
qualitative behavior of the field renormalization is already illustrated in 27J. It is essentially 
perturbative when k is larger than the mentioned threshold, with deviations of order 10~ 3 
from 1 (see Eq. ( 1351) ). Then, when k decreases, a bump appears for values of <fi around the 
separation point of the two regimes of V^. When k approaches zero this bump grows, dropping 
very rapidly to the standard perturbative value at the separation point. In Fig. [3] the field 
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Figure 4: Enlarged details of Fig. [3j Left panel: Zt (continuous) and Zl (dotted) for N = 4 at 
k = 0.02. Right panel: Z k (continuous) for N = 1 at k = 0.02. d = 4, A = 10 and A = 0.1 for all 
curves. 



renormalization in N = 1 and N = 4 (and c2 = 4) are shown. In the left panel of Fig. [3] we 
collected Z k (<f>) for N = 1, at k = 0.5 and k = 0.02 together with Z T for iV = 4, at k = 0.3 
and k = 0.02. In the right panel Zl for N = 4, is plotted at = 0.3 and k = 0.02 with a 
much larger scale on the y-axis because of the very high peak of Zl. 

Before going on, it must be remarked that the peak of -Zfc(0) for TV = 1 in Fig. [3] is much 
smaller than that observed in 27| (whose value was around 10), where the proper time RG 
flow with sharp cut-off on the proper time variable was used. That particular version of flow 
equation, due to its exponential nature, although very rapidly converging and very accurate 
in computing the critical exponents of the Ising universality class, was criticized in 25j for 
not well reproducing the details of the singularity in the second derivative of the potential in 
the broken phase. According to this criticism, the large size of the peak found in 27| could 
be addressed to that version of the flow equation which artificially amplifies the behavior of 
Zk close to 4>. On the other hand the approach followed in this paper is expected to be more 
accurate at quantitative level and therefore more suitable to describe this particular feature. 

Going back to the curves of the field renormalization, beside the size of the peaks other 
differences are illustrated in Fig. |4j where an enlargement of the lowest part of Fig. [3] (only 
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Figure 5: Zk(<p = 0) for N = 1 (black circles •) and Z £,((/> = 0) = Zt(4> = 0) for N = 4 (crosses X ) 
and Zi((j)) for TV = 4 ( plus +) at various values of k. In all cases d = 4, A = 10 and A = 0.1. Fits 
to these data are also plotted (see text). 



those curves at k — 0.02 are plotted) is shown. We note that both of the N — 1 theory 
and the transverse component Zt of the N = 4 theory show a very sharp drop followed by a 
drastic change of slope. However, by comparing the scale of the x-axis in the left and right 
panel of Fig. HI it is evident that Zy, with N = 1 is much steeper than Zt with N = 4 and, in 
particular, the drop observed in the N — 1 case is very close (within the numerical accuracy 
employed) to a discontinuous jump. Instead, the longitudinal component Zl of the N = 4 
theory smoothly approaches the perturbative regime for <fi ~ <fi and this, as shown below, 
yields a large (if compared with the perturbative correction), finite prediction of Zl(4>)- 

The structure of the flow equation of the field derivative is more complicated than the 
equation of the potential and even an approximate analytic solution is missing. However, 
from our numerical investigation we are now able to get some hints on the dependence of Z 
on various parameters such as the running scale k or the bare coupling A or the UV cut-off A. 
We concentrate on those values of <ft where Z can be reasonably extracted, and therefore we 
shall examine the origin = and the point 0. We avoid instead to analyze the peak of Z 
which is often affected by large fluctuations. 

We start by observing that in our computation in Fig. HJ for N — 1 the minimum is at 
4> = 1.285 and one finds Zk ~ 1.003 for (f> > 1.284 at k = 0.02, so that a typical perturbative 
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Figure 6: Z\~{^> = 0) for N = 1 (black circles •) and Z^{(j) = 0) = Zx(4> = 0) for N = 4 (black squares 
■) vs. y = 100 A at k = 0.02 and A = 10. Z k ((f> = 0) for N = 1 (stars *) and Z L (<p = 0) = Z T {(j) = 0) 
for N = 4 (crosses X ) vs. y = A at k = 0.02 and A = 0.1. All curves are for d = 4. Fits to these 
data are also plotted (see text). 



value is observed at 0. Then one should conclude that Zk=o(<fr) ~ 1.003, unless at k = the 
sharp drop of Z k hits 0, in which case no precise determination of Zk{<p) would be possible 
within our numerical accuracy. The data collected when k approaches zero do not exclude 
either of these two alternatives. In the case of the N = 4 theory one has = 0.819 and 
Zj< ~ 1.003 for > 0.816 at k — 0.02 similar to what happens in the N = 1 case with the 
difference that for N = 1 the jump could reasonably turn into a discontinuity, as already 
noticed, while in N = 4 the less steep drop reported in the left panel of Fig. |4] together with 
the smooth behavior of the potential around observed in Fig. |2j suggest that Z? could not 
develop a discontinuous gap at k = 0. On the other hand, as it is evident from the left panel 
of Fig. HJ the renormalization of the longitudinal field is more regular and the large peak 
decreases so smoothly that at one finds Zl = 1.324 at k — 0.02. Then, in this case it is 
possible to study the evolution of Zl{4>) as a function of k and eventually extrapolate its value 
at k = 0. 

In Fig. |5] we show some values of Z k ((f> = 0), Z L (<p = 0) = Z T ((p = 0), and of Z L ((j)), 
obtained at small values of k. In all cases A = 0.1, A = 10 and d = 4. For the points at = 
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Figure 7: Z L {4>) for for N = 4 (black circles •) vs. y = 100 A at k = 0.02 and A = 10. Z L (4>) for 
N = 4 (crosses X ) vs. y = A at A; = 0.02 and A = 0.1. In all cases d = 4. A fit to the latter set of 
data is also plotted (see text). 



we found that the functional form 

f(k) = ak 2 + b (36) 

provides excellent fits to the data with a = —0.079, b = 1.101 for N = 1 and a = —0.216, 
b = 1.094 for N = 4 and the corresponding plots are also shown in Fig. [51 As in the case of 
the potential, we observe here a quadratic dependence on the running scale k. 

Let us now consider Zl(4>). In this case the two expressions 

0.0014 



Z L = 1.0056 + — — and Z L = 1.0131 + 0.4225 e ~ liMbbk (37) 

k 2 + 0.0040 V ; 

fit equally well the data and we report both curves in Fig. [5j Their extrapolation at k = 
yield two different (of about 7%) values of Zl- However, remarkably, the fitted points do not 
suggest a divergent behavior of Zl(4>) at k — 0, but, on the contrary, indicate finite values not 
very distant from that found at k = 0.02 . We take this as a general indication that for this 
variable no strong modifications should be expected in the range < k < 0.02. 

In Figure [6] we show Z k (<fi = 0) and Z L ((f) = 0) = Z T (cf) = 0) versus y = 100 A , at k = 0.02 
and A = 10; we also show Zk(4> — 0) and Zl(<P = 0) = Zt(4> = 0) versus y — A at k — 0.02 
and A = 0.1. Except for the black circles, the other points show a rapid bend downward 
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for large y and this is due to the approaching of the critical value of A or A, which signals 
the transition to the symmetric phase. Apart from these points the two curves associated to 
circles and squares have the form g(X) = aX 2 + bX + c which is a typical expansion in powers 
of the coupling (for completeness we find a = —1.4823, b = 0.3855, c = 1.0700 for the squares 
and a = —1.4245, b = 0.5214, c = 1.0626 for the circles), while the points related to the field 
renormalization dependence on A display a linear behavior in a large range of the UV cut-off. 

Finally in Fig. we show two plots of Zl(4>), one versus y = A with A = 0.1 and the 
other versus y = 100 A with A = 10. All points are collected at k = 0.02. We observe a small 
dependence on the coupling with a small increase around A = 0.1 followed by a drop when 
approaching the critical line. Instead, we find a much stronger dependence on the cut-off with 
the field renormalization that reaches 3.5 for A = 3. By excluding the two points at large A, 
which are closer to the critical line, the other points are in good agreement with the curve 
Zl = 3.694 — 0.0562 e°' 401A , which can obviously be expanded in a polynomial. 



IV. CONCLUSIONS 



In this paper, we used the functional renormalization group to analyze the renormalization 
function of the scalar field in the N — 1 theory and of the longitudinal and transverse compo- 
nents in the 0(4)-symmetric theory in four Euclidean dimensions in the ordered phase where, 
due to spontaneous symmetry breaking, a non- vanishing vacuum expectation value <p ^ is 
generated. In particular, the approximation scheme of |8l-ll0| on the RG flow equations with 
the cut-off in Eq. ( 1261) provides a set of numerically stable coupled PDEs for V and Z, which 
allow us to approach the physical limit k = and to determine the main features of the field 
renormalizat ion . 

We found that a large non-perturbative enhancement of Z occurs in the range of <fi cor- 
responding to the observed flattening of the potential. This effect in Z had already been 
observed in quantum mechanics 24J and its large value in the classically forbidden region has 
also a direct interpretation in terms of suppression factor of the tunneling probability of a 
wave packet between two vacua ITJ, |47| . However, in field theory the numerical resolution of 
the equation is more problematical and a first indication of a large field renormalization for 



the N = 1 case was obtained in 27j. Our analysis qualitatively confirms those findings, but 
the maximum value of the peak found here is largely reduced with respect to 27 1. 

In particular we found that for N = 1 the potential Vj~ develops a non-analytic behavior at 
4> in the limit k — > and Zj. presents a peak followed by a sharp drop that occurs definitely for 
<p < at k = 0.02 (which is the typical lower limit reachable by our numerical routine). If this 
should persist even at k — 0, then Zk(<f)) would be substantially perturbative; otherwise if at 
k = the drop should occur at <p, then we were not able to obtain a reasonable determination 
of Zk(4>). The data collected close to k ~ do not exclude either of these two alternatives. 
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In the N = 4 theory the numerical analysis indicates that the potential does not develop 
a discontinuity in its second derivative with respect to <fr, as in the N — 1 case, and even the 
transverse component Zt shows a less sharp plot than for N = 1, suggesting in this case a 
continuous behavior at <fi and k ~ 0. At the same time the longitudinal component Zl has a 
very high peak but the curve is so well behaved that it is possible to extract Zl(4>) at various 
k and its extrapolation down to k = remarkably leads to finite values. 

The dependence of Z on the UV cut-off and on the bare coupling has also been investi- 
gated and this provides an indication on the typical numerical range spanned by the field 
renormalization, although it must be recalled that, quantitatively, these results are affected 
by large uncertainties because even a very small numerical error in the estimate of <f) induce 
large modifications in the determination of Zl. 
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